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Abstract. 

We present the first observational confirmation for an extended circumstellar dust disk around e Eri. The observations were 
obtained with the bolometer array SIMBA at the 15 m radio telescope SEST in Chile and measure the dust continuum at 
1.2 mm. The emission, with a total flux of 21.4 mjy and a rms of 2.2 mjy/beam, is resolved to a deconvolved size of 27'.'4 
which corresponds to 88 AU. No clear indication for a ring-like disk structure is seen, possibly also due to the telescope's large 
beam width of 24". Models of the object's spectral energy distribution from IR to mm-wavelengths show that the emission can 
also be explained by a simple disk model. We further demonstrate the strong influence of noise and propose to be cautious with 
interpretations of the ring substructure. 
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1. Introduction 

Circumstellar disks are an ubiquitous product of the forma- 
tion of stars and planetary systems. In these disks planets are 
believed to form from planetesimals that in turn are created 
through coagulation of dust and gas (Beckwith et al. 2000). The 
final stages of this process remain to be characterised but ob- 
servations suggest that previous protoplanetary accretion disks 
turn into debris disks at the end of planet formation. 

Photometric measurements with the IRAS satellite (Aumann 
I1988 1 provided first hints for dust around the ~800 Myr old star 
e Eri, a K2V dwarf at 3.22 pc distance to the Sun. From varying 
sub-mm flux densities in single-beam observations with differ- 
ent beam sizes (Chini et al. 1990 Chini et al. 119911 Zuckerman 
& Becklin 1993 ), an extended - possibly ring-like - disk struc- 
ture was assumed. In SCUBA observations at 850 /vm Greaves 
et al. ( 1998 ) found the dust emission peaking in a ring-like dis- 
tribution between 35-75 AU distance from the star. The cen- 
tral cleared region and an asymmetrical substructure in the ring 
were explained by the influence of orbiting planets. Now a 
planet with a semi-major half axis of 3.4 AU is known (Hatzes 
et al. 120001 1. which for more than a decade was controversial 
because of the star's high chromospheric activity. Quillen & 



Thorndike ( 2002 1 explained the ring substructure with a further 
hypothetical planet at 40 AU distance. CO line observations at 
the Swedish-ESO Submillimeter Telescope (SEST) concluded 
that the disk of e Eri is likely devoid of any gas (Liseau [T999> . 



With ISO no features from dust were found in mid-IR spectra 
between 6 and 12 //m, and in maps at 60 and 90 fj.m no resolved 
emission was seen (Walker & Heinrichsen 2000 1. A partially 
resolved region in IRAS data, quoted by some authors, was 
thus not confirmed. Several attempts to reveal the disk (in scat- 
tered light) with coronographic imaging at near-IR and optical 
wavelengths have been tried without success, including HST 
observations. Recently, Li et al. (|2D03 ) calculated a disk model 
from the IR to sub-mm spectral energy distribution (SED) of 
e Eri, but different disk and grain parameters were used than 
in our simulations. Furthermore, our calculations now also in- 
clude observations at 1 .2 mm. An overview of the fluxes used 
in this work is given in Tabled 

We present the first observational confirmation of an extended 
dust disk around e Eri found with the Sest IMaging Bolometer 
Array (SIMBA). With results from numerical simulations for 



Table 1. Fluxes of e Eri used in this work. The mid-IR IRAS 
fluxes contain also contributions from the stellar continuum. In 
the model SEDs in Fig.^this is separated from the dust emis- 
sion. 
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* Based on observations collected at the European Southern Obser- 
vatory, La Silla, Chile (71.C-0001) 



Wavelength [/mi] 


Flux [mJy] 


Error 


Reference 


12 


9520 


4% 


IRAS 


25 


2650 


6% 


IRAS 


60 


1660 


8% 


IRAS 


100 


1890 


9% 


IRAS 


450 


185 


103 mJy 


Greaves (1998) 


850 


40 


3 mJy 


Greaves (1998) 


1200 


21.4 


5.1 mJy 


this work 
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the target's SED we set constraints on dust properties and dis- 
cuss whether a continuous disk with an inner radius given by 
the dust sublimation temperature could be considered alterna- 
tively. 



2. Observation 

The observations of e Eri were obtained in July 2003 with the 
37-channel bolometer array SIMBA at the SEST in La Silla 
(Chile) using the fast scanning mode at 250 GHz (A = 1 .2 mm). 
Skydips were taken every 2 to 3 hours to correct for the at- 
mospheric zenith opacity, which was between to ~ 0.23 and 
to ~ 0.37. We followed the SEST pointing model (rms ~ 2'.'5) 
and checked the telescope pointing regularly. No significant de- 
viation from the pointing model was seen. The flux calibration 
is based on observations of Uranus. In total, 180 maps of e Eri 
were integrated which correspond to a total on-source time of 
~15 hours. 

Data reduction was performed with MOPSI 1 and included 
despiking, baseline fitting, suppression of the correlated sky 
noise, opacity and gain-elevation correction as well as co- 
adding the 180 single maps to a final one, which is shown in 
Fig-El The residual root mean square error (rms) amounts to 
2.2 mJy/beam. 

Observing at 1.2 mm with the 15 m dish of the SEST results in 
a nominal beam size of 24". This value can slightly vary with 
the elevation and other authors occasionally report 23" or 25". 
In principle, the beam size can be derived from a point source 
like, e.g., the planet Uranus. Averaging over various telescope 
positions we obtained a beam width of 24"2 + 072. The re- 
solved emission of e Eri can be fitted by an extended, circular 
Gaussian intensity distribution with a full width half maximum 
(FWHM) of 36"4. This fit and its centering is achieved with a 
routine in MOPSI. Relative to the beam width we resolve the 
dust emission around e Eri consistent with a deconvolved total 
extent of 27'.'4 corresponding to 88 AU. An elliptical shape of 
the intensity distribution - indicating an inclination of the disk 
- cannot be excluded, but is not apparent from our measure- 
ments. Greaves et al. ( 1998 1 reported i = 25°. At 1.2 mm, after 
subtracting a photosphere contribution of 0.9 mJy, the total disk 
flux amounts to 2 1 .4 mJy with a peak flux of 9.7 mJy/beam and 
2.2 mJy/beam noise. As can be seen in Fig. |2] (central panel), 
the peak of the Gaussian lies 7" southwest of the star. This is 
not necessarily a physical effect, since the accuracy of the cen- 
tering routine is within 5" (corresponding to 2.5 pixel in this 
figure). 



3. Numerical simulations 

The disk model is derived from the target's SED with the fluxes 
shown in Tabled Since especially the mid-IR fluxes contain a 

1 MOPSI has been developed and is maintained by R. Zylka, 
IRAM, Grenoble, France 



SED of Eps Eri 




A [fim] 

Fig. 1. Modelled SED of e Eri. The simulations were calculated 
with a grain size distribution and varying minimum grain size: 
flmin = 4//m (dashed line), a m ; n = 16 pan (dash-dotted), a m i n 
8 fim (solid curve, best fit). The stellar photosphere is shown as 
a dotted line. 

contribution from the stellar continuum, this is treated sepa- 
rately in the simulated SEDs. 

We apply a model of an optically thin disk with the simple 
radial density profile n{r) ~ r _1 (for the calculation of the 
dust temperature distribution see Wolf & Hillenbrand 2003 
Eq. 5). The following stellar parameters are used for e Eri: 
T* = 4700 K, R+ = 0.8 R Q and L* = 0.33 L (e.g. Soderblom 
& Dappen [l989> . We select a circular disk with an extension 
from Ria ~ 0.025 AU (corresponding to the radius given by the 
dust sublimation temperature) up to the observed R oM = 88 AU. 
Other choices for the inner disk radius and their consequences 
are discussed in the next chapter. 

Spherical homogeneous dust grains were used with radii a ac- 
cording to a size distribution n(a) ~ cT 3 - 5 . This describes the 
equilibrium size distribution resulting from a collisional cas- 
cade of dust (Tanaka et al. 1996 1 and has often been used 
for the interstellar medium (ISM) and modelling of circum- 
stellar disks. The best approximation to the observed SED is 
obtained for a minimum grain size a m [ n = 8 fim which is dis- 
played in Fig.[0 SEDs for a m i n = 4 /mi and a m i n = 16/im are 
shown for comparison. We use an upper grain size of 1 mm. 
The disk mass, given in the next section, refers to this up- 
per size limit, since we have no direct observational hints for 
cm-sized grains. Simulations are performed with a dust mix- 
ture deduced from the ISM: 62.5% silicate grains (accord- 
ing to Draine & Lee 119841 and Weingartner & Draine 120011 
mixed with 37.5% graphite. The dust grain density amounts to 
p(silicate) = 2.7 g crrT 3 and p(graphite) = 2.24 g cm -3 . 

4. Discussion 

Disk temperature and dust mass: 

Submillimeter observations probe the cold dust towards the 
outer part of the disk. From a fit to the observed fluxes be- 
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Fig. 2. Similar to the SCUBA observation our image sampling is 2". We find a central cavity and the emission peaking in certain 
spots along a ring-like distribution (left panel). However, subtracting a circular Gaussian intensity distribution with a FWHM 
of 36'.'4 (central panel) only leaves residues of pure noise (right panel). Thus we conclude that the substructure in our map was 
likely caused by remnant noise. The star marks the location of e Eri, corrected for its proper motion (-976.36 x 17.98 masec/year, 
cf. the SIMBAD database or Perryman et al. 119971 . In the left and central panel the intensity is scaled logarithmically from 5 to 
10 mJy, solid contours represent regions with an accuracy of 3 resp. 4 cr (rms = 2.2 mJy/beam), while the dashed contour (which 
nearly conincides with the outer margin of the greyscaled region) corresponds to 2.5 cr. The right panel shows -2 cr resp. +2 cr as 
dashed contours and +1.5 cr are symbolised by dotted lines. 



tween far-IR and mm-wavelenghts we determine for the cold 
dust component a temperature of 50 K using (3 — and 
k{\.2 mm) = 1.7 cm 2 g _1 (a value often used for Vega-type 
disks, e.g. Sylvester et al. 2001 1. With the numerical simula- 
tions described above we obtain a dust mass of 5.2 x 1O~ 9 M 
= 1.7 x 10~ 3 M ffi , which lies about a factor of three under the 
lower mass limit given in Greaves et al. J1998l >. 

Dust ring or a continuous disk: 

We are able to explain the observed SED with a simple disk 
model, an inner disk radius given by the dust sublimation tem- 
perature and a minimum of free parameters, i.e. without fur- 
ther assumptions. However, we do not exclude other disk mod- 
els. With the limited resolution of the SEST beam only a very 
sharp ring would have been detectable. Even for SCUBA the 
reported central cavity is close to the resolution limit. To ex- 
plain however the mid-IR IRAS fluxes, some amount of hot 
circumstellar dust must be present at regions relatively close to 
the star. From modelling the observed SED no statement can 
be made to which extent the inner disk region may be depop- 
ulated of dust. Due to the large parameter space in any dust 
disk model (disk density profile, dust grain radii, grain size dis- 
tribution, dust mineralogy, grain shapes, etc.) several solutions 
may exist to interpret the observed SED. Li et al. (2003) find 
their model in good agreement with Greaves et al. ( 1998 1 for 
disk radii > 28 AU, while an additional "zodiacal" dust compo- 
nent would be required to explain the possible emission from 
the inner parts of the disk. This inner region was not clearly 
resolved by previous mid- and far-IR satellite cameras, but can 
be better studied with upcoming SIRTF observations (see Li 
et al. 120031 for a discussion on this). An alternative approach 
would be interferometric observations with the IRAM Plateau 



de Bure Interferometer operating at 1.3 mm wavelength with a 
resolution up to 1". 

Influence of noise on the observed disk substructure: 

Fig-E](left panel) shows our final co-added map of e Eri. The 
pixel sampling with 2" is justified because of SIMBA's im- 
age sampling of 15 ms in azimuthal scan direction. With 80"/s 
scanning speed this results in a spatial sampling of 172. In right 
ascension and declination we obtain a similar value due to the 
large coverage of maps taken at different phase angles. Our fi- 
nal image resembles a disk with a cleared central region and 
substructures in the ring. However, we can show that these de- 
tails are most likely caused by remnant noise. The resolved 
emission of e Eri can be fitted by an extended Gaussian in- 
tensity distribution with a FWHM of 3674 (central panel). A 
subtraction of this fit from the final map only results in pure 
residual noise, while a real central cavity should cause a clear 
negative peak after subtraction. The occurence of a minimum 
exactly at the position of the star leaves speculations whether 
this is a contribution of a real gap. However, this minimum is 
even less prominent than the noise minima above and below 
the star's position. Similarly, the knots in the outer disk part 
can be explained by noise maxima. To verify this, we simu- 
lated a ring-like disk with artificial noise but without intrinsic 
ring substructures, and obtained a comparable result. 

The knots in the SCUBA ring were attributed by some authors 
to be traps of dust due to orbital resonances of Jovian planets. 
However, most of these are smaller than the resolution of the 
SCUBA beam. We speculate that these structures may be ex- 
plained similar to our observations. 
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5. Conclusions 

Since several publications calculate planetary parameters from 
the knots in the e Eri ring structure, we intended to perform 
the first confirmation of this disk shape. Our observation with 
the bolometer array SIMBA operating at 1.2 mm is the first 
verification of an extended disk around this star. However, these 
data do not provide clear indications for the previously reported 
ring substructure. In our case the knots can be explained with 
remnant noise effects. We are able to model the observed SED 
with a simple, gap-free disk, but do not rule out the existence 
of a possible central cavity. 
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